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CN ■ QCD predicts the existence of a new partonic Yukawa phase above a critical temperature 

Oh' Tc ~ 150 MeV. I review some of the key observables which will be used in heavy ion 

reactions to search for this new phase at RHIC. Systematics of collective observables 

{Nch,Ej_), flow patterns, meson interferometry, jet quenching, and J/ip suppression are 






^ I discussed. This talk is updated with first data from RHIC 
^ . 1. Yukawa in QCD 

cn ■ 

C^^ . In 1935 H. Yukawa|l| proposed a theory of nuclear forces based on the exchange of a 

^ ! massive boson 

O 

§. Vy{r,m) = a^ff^ (1) 

He estimated that m ~ 100 MeV to account for the short range ~ 2 fm of the nuclear 

force, and in 1947 Powell confirmed his theory. Yukawa's meson theory forms the basis 
^ I for the current effective theory of nuclear forces]^. 
C ■ In the field of electrolytes, Debye and HiickleQ] had already discovered the Yukawa 

potential in a completely different context. The polarizability of the medium in the 
^ ', presence of an external charge density leads to a non-linear self consistent equation with 
H [ the same solution as eq. (|l]), but in that case the effective mass is the Debye electric 

screening mass 

^^' = AnJ2q'\n,\/T . (2) 

A conductive medium transforms Coulomb power law forces into the Yukawa form. 

In nuclear theory, the Yukawa meson mass results from the finite gap of the elementary 
excitations (pions, ...) of the physical QCD vacuum. In this talk, I discuss current efforts 
to try to drill a perturbative hole into the nonperturbative vacuum using RHIC to see 
the breakdown of Yukawa's meson theory. As we review below, in the deconfined, chirally 
symmetric phase of QCD at high temperatures, the Debye- Huckle mechanism transforms 
the mesonic Yukawa potentials into a color-electric screened gluonic Yukawa potential in 
a quark-gluon plasma (QGP). 

In QCD, the color potential between partons is approximately Coulombic at small 
distances due to the asymptotic freedom property of non-Abelian gauge theories. However, 
below a critical temperature, T^ ~ 150 MeV, the effective potential between the colored 
partons has a long range linear confining term Kr with a huge "string" tension, k ~ 1 



GeV/fm. In this confining phase of QCD, the heavy qq potential is well parameterized 
by the Liischer forni|^ 

VL(r,0) = -— + /£r , (3) 

r 

(as long as dynamical quark pair production is ignored). The Coulombic part, with 
strength a^ = tt/12, arises from the zero point quantum fluctuations of the string. As 
the temperature increases, but remains below the deconfinement transition, T < Tc, the 
enhanced fluctuations due to thermal agitation of the string modifies the effective potential 
into the approximate Gao form [^ 



VGir,T) = -^ 
r 



2 
1 tan"^(2rT) 



K i 1 tan ' 



3 ' n '2rT' 



r + ■ ■ ■ (4) 



The decrease of the effective string tension, k{T), predicted above has been measured via 
lattice QCD calculations P as shown in Fig.(|ID. However, the "measured" string tension 
is found to decreases faster than predicted in eq.(^) near the critical temperature. Note 
that in Fig.(|I|) T ^ T,exp(ll/6(/? - /?J) with p^ = 4.0729 for this lattice. 

For temperatures above Tc, the lattice data in Fig.|l] reveal the predicted Yukawa phase 
of QCD is predicted. The heavy quark potential mutates into a short range generalized 
Yukawa form, which on the lattice is fit with the form 

t'anr,d) = -i^e-<-W^. (6) 

Note that di = 2,/xl(T) = 2mE(T) correspond to the ideal Yukawa form. The perturba- 
tive thermal QCD chromo-electric Debye mass niE = /x(T)/2 is 

™.(T)^.(T)T(f.^)"'^ (6) 

for Nc colors and Np flavors. For Nc = 2,Nf = in Fig. (|l]), we expect ^l{T) = 
2mE = 1.6g{T)T as shown by the solid line in Fig.(|^). The fits to the lattice QCD data 
in Fig.(|I| b) from show that for T > 2Tc fii ~ 2.5T is not far from the pQCD estimate. 
However, the exponent di ^ 1.5 is significantly below the value 2 expected from pQCD. 
An even more striking nonperturbative deviation is seen near Tc, at which point d < 1 
and fi ~ T/2. This suggests a rather long range interaction that may be the precursor of 
the confinement transition. 

The predicted thermodynamic properties of the QCD Yukawa phase are shown in Fig.^ 
from ref. 1^ for 2 flavor 12^ x 6 IQCD. A present caveat about all lattice results is that 
the pion is still too massive to allow contact with the "known" thermodynamic properties 
of ordinary hadronic /nuclear matter below Tc. Two striking features of Fig.(|]) suggest 
two key observable signatures of this phase transition in nuclear collisions. First, the 
entropy density clT) = (e + p)/T increase very rapidly with T in a narrow interval 
AT/Tc < 0.1. Second, the plasma becomes extremely soft p/e <^ 1 and c^ -C 1 near T^. As 
we review below, the first feature may lead to time delay (the QGP stall) measurable via 
hadronic interferometry. The second feature may lead to interesting non-linear collective 
flow observables in nuclear collisions. The experimental verification of these fundamental 
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Figure 1. The heavy quark potential in the confined phase of SU(2) quenched QCD (top) 
compared to that in the deconfined Yukawa phase above Tc (bottom). Results are for a 
32^ X 4 lattice from Karsch et al0. Eq.(^) fits the confined lattice "data" well, but the 
QCD string tension decreases more rapidly near T^. For T > T^ the potential is screened by 
the deconfined gluons (in this quenched calculation) and acquires the generalized Yukawa 
form (D. Here T ^ Tcexp(ll/6(/3 - /5J) increases with (3 and (3^ = 4.0729. 




Figure 2. The chroino-electric Debye screening mass, fi{T) = 2mE{T) and the effective 
exponent, d{T), of the effective Yukawa potential in the deconfined phase versus T/Tc is 
shown from Karsch et al for SU(2) quenched QCD. The hues show expected dependence 
based on thermal pQCD. Note that near Tc, the range 2//i is much larger than predicted 
by pQCD and that (i ~ 1 implies an especially long range interaction there. 
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Figure 3. Thermodynamic energy density {e/T'^ top curves left), pressure {3p/T^ lower 
curves left) and speed of sound squared (right) from lattice QCD (2 flavor 12^ x 6) from 
the MILC collab 0. The curves are zero quark mass extrapolations. Note the rapid 
reduction of the pressure and speed of sound as Tc is approached from above. 



predictions of QCD is one of the primary goals of the heavy ion experimental program at 
Brookhaven and CERN. In the following sections, I review first how deep into the QGP 
phase RHIC may be able to reach, and then discuss several signatures that will be used 
to test the QCD predictions in such experiments. 

2. Initial conditions in A+A 

In order to see the partonic Yukawa phase, we must first heat the vacuum to about 100 
times the energy density of nuclear matter. At SPS energies the low p_j_ physics makes 
it impossible to predict the initial conditions. However, with increasing energy pQCD 
begins to provide increasingly more reliable theoretical basis for predicting those initial 
conditions. For highly boosted nuclei with Ecm > 100m j^j, time dilation effectively freezes 
out the quantum chromo fluctuations inside the nuclei while the two pass through each 
other. Au beams at collider energies can be thought of as well collimated, ultra dense 
beams of partons. This ( chromo Weizsacker-Williams0) gluon cloud contains a very 
large number, Ga{x,Pq) ~ A/x^^^, of almost on-shell collinear gluons with longitudinal 
momentum fraction x = po/Ecm -C 1. As the clouds pass through each other, many of 
the (virtual) partons scatter and form a dominantly gluonic plasma on a very fast time 
scale l/p± -C 1 fm/c. The number of gluons pairs (mini-jets) extracted from the nuclei 
by this mechanism at rapidities t/i and transverse momentum ±k_j_ can be calculated in 
pQCD from the well known expression ||T0|-[T3[ 



f^f^'S = Kx,Ga{x„ kl)x,Gs{x2, kD^^^^TMb) , (7) 

dyidy2dk\^ akj^ 

where Xi = x±{exp{yi) + exp{y2)) and X2 = x±{exp{-yi) + exp{-y2)) , with x± = k±/y/s, 
and where the pQCD gg -^ gg cross section for scattering with t = —k'j_{l + exp(?/2 ~yi)) 
and y2 — yi = y is given by 



da39 _ 9 47ra2(l + e^ + e-^)^ 



dt 8 ki {ey/^ + e-y/^Y 



(8) 



For atomic numbers, A ^ 1, the geometrical amplification, Tab(6) ~ 30/mb, enhances 
by orders of magnitude the gluon density relative to pp. The factor i^ ~ 2 approximates 
next to leading order contributions. 

For symmetric systems, A + A, with G^ ~ AG, the inclusive gluon jet production cross 
section is obtained by integrating over ?/2 with yi = y and fcj_ fixed. To about 50% accu- 
racy, the single inclusive gluon rapidity density in central collisions is approximately |T^ 

dN A'^ dcT^d 

i;^.A^27V,(.,..).,G(...<)^o.AV» (9) 

where Ng(x±,t) = J^ dxG{x,t). This copious mini-jet production mechanism is believed 
to the dominant source of gluon plasma production at RHIC and higher energies. 

Recent upper bound estimates of the total gluon rapidity density as a function of the 
CM energy from EKRT[|l^ are shown in Fig.(^. The differential yields are integrated 



down to a transverse momentum scale pq ^ 1 — 2 GeV. This scale separates the "soft" 
nonperturbative beam jet fragmentation domain from the calculable perturbative one 
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Figure 4. Mini-jet initial rapidity density of gluons, A^j, produced in central A+A col- 
lisions as a function of the pQCD cutoff, po? from |jlO[|. The dashed "saturation" curve 
corresponds to (pq-R^). (A^j is dimensionless; [mfe] is a typo). The magnitude of possible 
hydrodynamic transverse energy loss due to longitudinal work in an ideal p = e/3 quark 
gluon plasma is shown by the difference between Eti and E't/ in the right panel. 



above. The curve marked saturation |T^ is an upper bound marking the point where the 
transverse gluon density of mini-jets becomes so high that the newly liberated gluons 
completely fill the nuclear area, i.e., dN/dy ^ PqR^- A current hypothesis is that at 
that point, higher order gluon absorption may limit the further increase of the gluon 
number. At RHIC energies EKRT estimates yield up 1500 gluons per unit rapidity. Our 
conservative estimates together with X.N. Wang |lI^JT5|l only gives < 500 gluons per unit 
rapidity when initial and final state radiation is also taken into account. This lower 
number is obtained with a fixed po = 2 GeV, that we found to be necessary in order 
to reproduce pp, pp and lower energy BA data using the HIJING event generator 



The initial energy density reached in such collisions can be estimated using the Bjorken 
formula 



<ro) 



dEj 



(10) 



ttR!^tq dy 

For Po ~ 1 - 2 GeV, dEr/dy ~ 400 - 2000 GeV, and therefore e(ro ~ 0.5 fm/c) > 10 
GeV/fm^ should be easily reached at RHIC, well inside the the deconfinement phase of 
QCD. At SPS energies, on the other hand, our estimates indicate that nuclear collisions 
may just reach the transition region and cool below it in a very short time. 

Recently CERN issued a press release [^ claiming that " We now have evidence of a new 
state of matter where quarks and gluons are not confined." As discussed below, I disagree 
with the claim that deconfinement has been observed. What I find compelling is that some 
form of matter, much denser than ever studied before, was created. Inferences about the 
role of quark and gluon degrees of freedom cannot be drawn because at the relatively low 
momentum scales accessible at SPS energies, the quarks and gluon degrees of freedom in 
the dense matter are mostly not resolvable. Even at the highest transverse momentum 
measured, the pion spectra were shown to be very sensitive to nonperturbative model 



assumptions regarding the role of intrinsic momenta and soft initial state interactions! 17 



The dynamics of the non-perturbative beam jet fragmentation and hadronic final state 



interactions simply cannot be disentangled at SPS energies[lS]. While there are many 
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Figure 5. Charged particle rapidity density per participating baryon pair from [^ is 
shown versus the cm energy per baryon. The PHOBOS data[|19| (filled triangles) for the 
6% most central Au+Au are compared to pp and pp data (open symbols) and the NA49 
Pb + P6(central 5%) data H1J1NG1.35 (solid) and EKRT (dot-dashed) predictions are 
also shown. HlJlNG upper (lower) solid curves are without (with) jet quenching. The 
right panel shows the characteristic difference between the predicted dependence of the 
charged multiplicity on the number of participant baryons for HIJING (solid) and EKRT 
(dot-dashed). 



interesting signatures showing that dense matter was formed at the SPS (through the 
non-linear in dependence of several observables on multiplicity or A), the bottom line is 
that those data have said nothing about whether the QCD predictions in Figs 1-3 are 
correct or not. We need higher resolution, i.e. energy, to see the quarks and gluons in 
action in the new Yukawa phase. 



3. First RHIC Data on Multiplicity 

The first published data[|T^] from RHIC are shown in Fig.^. Amazingly, both HIJING 
and EKRT predictions are seen to be consistent with the central multiplicity density data 
in Au+Au at ^/s = 65, 130 AGeV as measured by PHOBOS [|^. However, as emphasized 
the scaling as a function of centrality will soon be able to differentiate between 



m 



20 



these two very different pQCD based estimates of RHIC initial conditions. In addition, the 
transverse energy systematics will provide an independent critical test of the hypothetical 
saturation picture |2^ as we discuss below. While it is obviously premature to draw 
any physics conclusion from these first data, the consistency of the increased activity 
per baryon predicted due to the onset of mini-jet activity beyond i/i > 100 AGeV is 



reasuring. A critical test of the pQCD framework used above will be the upcoming RHIC 
data at ^/s = 200 AGeV. HIJING predicts sensitivity to jet quenching at that energy. 

4. Barometric Measures of Collective Dynamics 

The simplest global barometer of collectivity in nuclear reactions is the A and energy 
dependence of the transverse energy and charged particle rapidity density. At RHIC 
energies, HIJING predicts that initial transverse energy density in central A + A collisions 
scales nonlinearly with A 

^.lGeVA-,l,,og^) (U) 

This leads to about 0.6 TeV per unit rapidity at the present ^/s = 130 AGeV energy. In 
EKRT[0, on the other hand, the gluon saturation hypothesis predicts an approximate 
linear A dependence of E± ~ A^^^, and a value several times that of HIJING. In Fig.(^), 
the initial gluon density however grows less than linear A^'^"^ in that model. The A^'^ 
scaling of HIJING with the fixed mini-jet scale po = 2 GeV scale. This results from the 
combined increase of the number of binary interactions as well as the increased fraction 
of the energy that originated from the mini-jets. The initial energy density in HIJING 
varies approximately as 

Co ^ 0.6 GeV/fm^ A^-^^ (1 + log ^) . (12) 

EKRT Q on the other hand predict e(r = l/psat) ~ 0.1/l°-^s°-3^ 

If local equilibrium is achieved and maintained, then hydrodynamics predicts that lon- 
gitudinal boost invariant expansion together with pdV work done by the plasma pushing 
matter down the beam pipe will cool the plasma and convert some its random transverse 
energy into collective longitudinal kinetic energy. For an equation of state, p = c^e, this 
cooling and expansion causes the energy density to decrease with proper time as 



e(r)=e(ro)(^)'^'^ (13) 

Entropy conservation leads to a conservation of rnij), where n is the proper parton 
density. At least if the matter is initially deep in the plasma phase, then (as seen in 
Fig. 3) longitudinal will be done with c^ k, 1/3. Consequently, the transverse energy per 
particle should decrease by a factor 2-3 before freeze-out as|21[] 

However, dissipative effects due to rapid expansion and finite mean free paths reduce 
considerably the effective pressurep2|. For the Bjorken expansion, the relaxation time, 
Tc = l/{aTn) oc r, then increases with time as n decreases. Numerical solution of 3+1 
D transport equations with pQCD cross sections, a^ ~ 2 mb, indicate that dissipation 
reduces the transverse energy loss for HIJING initial conditions rather significantly (see 
detailed comparisons in p3| , p4| ). 



One of the important experimental observations at the SPS is that at y/s ~ 20 AGeV, 
dE±/dy as well as dN/dy scale nearly linearly with wounded nucleon or participant num- 
ber (~ A^-^''') 1^. Simple Glauber wounded nucleon models reproduce very well the nearly 
linear correlation between Ej_ and the veto calorimeter (spectator) energy observed in all 
experiments at SPS|^. The physics implications of this scaling depends on what is as- 



sumed for the A dependence of the initial conditions. One view0] is that the initial e± 
scales in just the right way that after hydrodynamic expansion the final Et and e± always 
scale linearly with A. My view is that at the SPS the linear dependence arises from ad- 
ditive nature of soft beam fragmentation together with the absence of pdV work at early 
times. If the QGP transition region is just barely reached, as I believe, then the softness 
of the QCD equation of state with c^ -C 1/3 seen in Fig.(H) and dissipation conspire to 
prevent the dense matter from performing longitudinal work. However, it is impossible 
to tell from the data whether this observed null effect is due to a low pressure Hagedorn 
resonance gas of hadrons or to a low pressure lazy "plasma" with c^ -C 1. At RHIC en- 
ergies the plasma starts deep in the Cg ~ l/sqrt3 Yukawa regime, and longitudinal work 
should become observable in the Et systematics. 

5. Transverse Collective Flow 

An entirely different measure of barometric collectivity is afforded by the study of the 
triple differential distributions, dN/dyd'^'pj_. Already at sub- luminal Bevalac energies (< 1 
AGeV), azimuthally asymmetric collective directed and elliptic flow were discovered long 
ago. For non central collisions, 6 7^ 0, the asymmetric transverse coordinate profile of the 
reaction region leads to different gradients of the pressure as a function of the azimuthal 
angle relative to beam axis. This leads to a "bounce" off of projectile and target fragments 
in the reaction plane and to azimuthally asymmetric transverse momentum dependence of 
particles with short mean free at mid rapidities. This phenomenon has now been observed 
at both AGS and SPS energies as well[|^. It has been predicted to be there also at RHIC 



and LHC energies[2S]. In Fig.(H) the first STAR datap9| on the centrality dependence of 



the second Fourier component, V2, of the azimuthal flow patterns is shown: 

ocl + 2t;icos(0-0ij) + 2f2cos(2(0-0ij)) + --- . (15) 



dyd^p^ 



These data show that the azymuthal asymmtery is about twice as large at RHIC than at 
the SPS. Furthermore, the differential azymuthal flow, f2(p±), rises linearly with p_|_ up 
to 2 GeV/c. 

The important question is how this type of barometer could serve to help search for 
evidence of the QCD transition. In |^ the idea was proposed that one could use V2 
systematics to search for the predicted softening of the QCD equation of state. Hydrody- 
namics calculations lead to a factor of two smaller V2 for an equation of state with a soft 
critical point as in Fig 3 versus one in which the speed of sound remains 1/a/3. Searches 
for anomalous V2 dependence as well as Vi are underwayPT|. As with the global barome- 



ter, dissipation can of course also simulate a soft equation of state. In ref [^ we studied 
the dependence of f 2 on the transport parton cross section for RHIC conditions. The 
results are rather sensitive to the cross sections and initial densities as shown in Fig(B[). 
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Figure 6. (left) First STAR data ^^ on the centrality dependence of the azimuthal flow 
moment f 2 for Au + Au at 130 AGeV. Estimates of the corresponding spatial asymmetry 
are indicated by boxes. (Right) Predicted dependence of V2 at mid rapidity Au + Au 
at RHIC from parton cascade calculations ZPC[^] with HIJING initial conditions for 



agg = 3 and 10 mb. (Note V2 is the p± weighed form of ^2) 



For HIJING initial conditions, dissipation leads to a significant reduction of f 2 relative to 
hydrodynamics. 



6. The QGP Stall and Interferometry 

Hadron interferometry has been developed over the last several decades in heavy ion 
collisions into a precise tool to image the space-time region of the decoupling 4- volume. In 
Q it was proposed that a possible signature of the QGP phase transition would be a time 
delay associated with very slow hadronization. The plasma "burns" into hadronic ashes 
along deflagration front that moves very slowly if the entropy drop across the transition 
is large pi 



In |5§ 



we calculated the 2+lD evolution of a Bjorken cylinder with time 
and transverse coordinate using hydrodynamics with different equations of state. For an 
initial energy density ~ 20 GeV/fm^ a time delay up to a factor of two was found even 
for a continuous transition as long as the entropy jump was ~ 10. This QGP "stall" of 
the transverse expansion is due to the small speed of sound in the mixed phase. It should 
be readily observable in pion interferometry by comparing Rout and Rside radii. 

The time delay is a robust generic signature of a rapid cross over transition of the en- 
tropy density. In particular the "stall" is expected even for a smooth cross-over transition 
as long as the width AT/T^ < 0.1. However, its magnitude also depends strongly on 
the entropy drop across that region. Unfortunately, as noted before lattice QCD has not 
yet been able to resolved the hadronic world below T^ due to technical problems. If the 



11 



entropy jump in nature is less than a factor of three, then the stall signature would be 
much more difficult to observe. 

High statistics measurements of pion and kaon interferometry searches for time delay 
at AGS and SPS have come up empty handed thus far. No time delay has ever been 
observed in any nuclear reactions thus far. This could be due to (a) the absence of a 
large rapid entropy drop in real QCD, or (b) to unfavorable kinematic conditions at AGS 



and SPS energies. From the hydrodynamic calculations in |35], it was found that a large 
time delay signal requires a favorable initial condition with initial temperature ~ 2Tc. 
For too high an initial temperature (as at LHC) the large transverse collective flow that 
develops prevents a stall from happening. For too low a temperature (as at SPS) the 
time delay is suppressed because the system spends too little time in the mixed phase 
as a result of longitudinal expansion. Therefore, RHIC is the most likely energy regime 
where a QGP stall may be observable. Recent work |^ has also shown that high p± kaon 
interferometry is an especially sensitivity probe of time delay. If observed, a time delay 
signature would be one of the most powerful indicators of novel collective behavior that 
can only arise if the produced matter has an anomalous (softest point) equation of state. 



7. J/Psi Suppression 

In 1986, Matsui and Satz proposed an intriguing direct measure of the transmutation of 
the qq forces shown in Fig. 1. The idea was that J/ip can form in the vacuum because the 
confining Luscher potential can bind a cc pair into that vector meson. If that pair were 
placed in a hot medium in which the potential is screened, then above the temperature 
where the screening length is smaller than the J/ip radius, the cc would become unbound. 
The charm quarks would then emerge from the reaction region as an open charm DD 
pair. 

J/ip suppression was soon seen in 1987 in O + U reactions by NA38. Since then this 
smoking gun has (unfortunately) never stopped smoking! J/ip suppression seems to be 
as ubiquitous as the Yukawa potential. It was subsequently observed in p + A reactions 
as well. High mass Drell-Yan pairs, on the other hand, formed via qq -^ £i was observed 
to scale perfectly linearly with the number of binary collisions. This is because lepton 
pairs suffer no final state interactions and the quark initial state (Cronin) interactions are 
invisible in p± integrated DY yields. 

The "normal, conventional" J/^p suppression leading to a [AB)^'^ dependence, is pre- 
sumed to be due to a mechanism that is independent of QGP production. In Pb+Pb, 
NA50 observed an excess 25% suppression of J/tp above that normal expectation. NA50 
therefore claimed that this enhanced suppression is finally the sought after smoking gun 



of QGP format ion p7 



While the deviation from the empirical [AB)^'^ scaling is convincing, the dynamical 
origin of the effect is far from settled in my opinion. As shown by the several curves in 
Fig.d^), a ~50% drop in the J/ip yield as a function of Et is consistent with final state 
co-moving hadronic absorption. While the detailed wiggles are not reproduced, large 
theoretical uncertainties about several key dynamical ingredients preclude more precise 
comparisons at this time. It is important to emphasize that plasma scenario models suffer 
at least as large theoretical uncertainties as the so called "conventional hadronic" models. 
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Figure 7. "Evidence for deconfinement of quarks and gluons" claimed by NA50 in Ref . [^ . 
The curves on left show transport theory estimates for hadronic final state dissociation. 
The enhanced suppression relative to exp{—aNPoL) nuclear suppression observed in p + A 
and S + U is shown versus a very rough estimate of the unobserved initial energy density. 



Key uncertain elements include (1) the proper QCD treatment of cold nuclear absorption 
responsible for the nonplasma AB ' suppression observed in p + A, (2) the unknown 
hadronic M + ip -^ DDX reaction rates, and (3) the detailed density evolution, e{x,t). 

For example, only a schematic "octet model" of pre-hadronization cc exists at present 
to calculate the so called "normal" nuclear suppression. That this is model may not be 
suffiencient to adequately account for that component can be seen in the work of Ref. [^ . 
The observed J/tp production cross section (ignoring final state interactions) may be 
expressed as 



CAB^ipX 



daAB- 



tcc" cc—*ip 



{<f 



(16) 



where F is the formation probability of the ip from a cc that emerges from the cold nuclear 
target with an invariant mass q^ < 4M|,. If the pre-resonance cc pair multiply scatters in 
the nucleus, then q"^ —>■ q"^ + 5q^{apL) increases linearly with nuclear thickness as shown 
by the curve G in Fig.(P). This leads to an approximate exponential suppression that 
can account for the approximate Glauber nuclear absorption factor ansatz, exp(— cre//pl/). 
This Gaussian model can thus account for the observed {AB)^'^ scaling for light projectiles. 
However, it was shown in [^ that if the power law tails due to induced radiation in the 
medium are included (resulting from the multiple Rutherford rescattering of the color 
octet cc) , then an additional nonlinear suppression in the nuclear thickness L could 
result (curve P). This is because induced gluon radiation provides a way to increase the 
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Figure 8. Sensitivity of "pre- resonance" J/ip nuclear absorption to details of color octet 
models of cc interactions in a nuclear medium from Qiu et aliSH]. The curve marked P 



for power law accounts for the anomalous absorption seen in Pb + Pb and deviated from 
Gaussian Glauber-like G expectations. 



invariant mass of the pre-resonance cc that can further reduce the probability for the pair 
to fit inside the tp wavefunction. While this model is also rather schematic, I feel that 
it captures one of the possibly important elements of pre-formation transient dynamical 
effects in nuclei that may account for anomalous nonlinear suppression in the thickest 
nuclei. 

The claim of anomalous suppression cannot therefore rest merely on generic enhanced 
suppression in Pb. It rests on the possible existence of singular "step-like" structure of the 
suppression pattern. The evidence for "steps" is however the weakest link experimentally 
because a rigorous x^ test including the substantial systematic errors in the Et scale has 
yet to be performed. Nevertheless, if the "step-like" suppression pattern survives further 
experimental scrutiny, it would certainly be the most dramatic nonlinearity observed at 
SPS. Experimentally, the claims would carry considerable more weight if similar "step- 
wise" patterns were observed in other systems ,e.g. Xe + Xe suitably shifted in Et due to 
the expected smaller energy densities achieved. I would put such an experiment highest 
on the SPS priority list. At RHIC, the striking prediction by H. Satz at QM99[^, was 



that under RHIC conditions the higher energy density should lead to the same step-wise 
pattern in the lighter nuclear Cu + Cu interactions. PHENIX will provide a definitive 
test of this prediction in a few years. 

One significant inconsistency with the present plasma scenario interpretation is its 
failure to account for the observed Et dependence of the J/ip p± spectra in Fig.(^. 
Standard Glauber multiple collisions lead to a random walk in transverse momentum 
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Figure 9. The mean transverse momentum of the J/Tp |]5D| shows clear evidence of multiple 
scattering in the nuclear target as expected ||41|| from Glauber theory. Even in Ph + Ph 
the increase of the pt is understood from the same mechanism in contradiction 1 5^ to 
predictions based on the plasma scenario, where at high Et-, the surviving if) are expected 
to have been produced only near the nuclear surface regions |^|. 



that are expected enhance the J/^ transverse momentum as pi 
L 



{p±)ab 



ipi) 



PP + j^Pa 



(17) 



This is found to hold[^ experimentally in all reaction including Ph + Ph. In contrast, 
in the plasma scenario |Q, only those ip are expected to be observed which are produced 
near the surface, where the effective nuclear depth L is small. Thus, the prediction as 
shown in Fig(B|) was that the (p^) should begin to DECREASE with increasing Et- This 



was not observed [39, 40 1 . 



8. The High p±^ Frontier 

One of the most exiting new physics areas that will become accessible at RHIC is 
the study of high transverse momentum (short wavelength) jet probes. The rates of 
jet production and its fragmentation in the vacuum are well understood. The new jet 
physics at RHIC will be the study of partonic interactions at extreme densities through 
the phenomenon of jet quenching [^] . Final state interactions of a jet in a dense QGP are 
expected to induce a large radiative energy loss[^. In fact, BDMPS[^ discovered that 
non-Abelian energy loss is in fact non-linear as a function of the thickness of the medium. 
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Figure 10. Jet quenching at SPS vs RHIC from ref. |T^ compared to WA98 data. At 
RHIC the power law tail extends finally far enough above the nonperturbative "noise" to 
make jet quenching observable (long dashed in right panel). At SPS the steep high p±_ 
tails are too sensitive multiple soft scatterings. 



Tests of this and other aspects of non-Abelian multiple collision dynamical phenomena 
will soon be possible at RHIC^ . 

At lower SPS energies, this physics is out of reach because nonperturbative multiparticle 
production physics (e.g., soft multiple p± kicks) dominates as shown in Fig.(pi]|) from ref. 
]17| . HIJING happens to fit the WA98 data with or without jet quenching at SPS energies. 
While at the SPS no clean separation of soft and hard dynamics is kinematically possible, 
at RHIC energies, the high p_j_ power law tails of the single inclusive distributions stick out 
far enough above the soft physics "noise" that sensitivity to the form of the non-Abelian 
dE/dx is expected as shown in Fig.(|10|). This problem is closely related also to the problem 
of pre-resonance J/ip absorption discussed previously |3^. Understanding jet quenching 
is also necessary to develop more accurate covariant parton transport theories p3|,p4|. 



9. Summary 

The next Yukawa phase of QCD is awaiting discovery. The SPS data have provided 
many intriguing indirect hints that new physics is already operating near the transition 
region. Many claims and counter claims remain at SPS energies because both hadronic and 
partonic models have partial overlapping domains of validity. This "duality" is analogous 
to the problem of interpreting the R factor in e^e~ collisions below ^Js < 10 GeV. The 
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ratio of hadronic to leptonic production cross sections only reaches the magic 11/3 of 
pQCD above the threshold regions for ss, cc, bb. Near the threshold regions conventional 
vector mesons models of nonperturbative hadronic physics provide an adequate "dual" 
description of the physics. The simpler continuum partonic description is only applicable 
far above that region. Similarly, reactions at SPS energies are near the threshold region 
where resonances become too broad due to multiple interactions. In nuclear reactions, the 
continuum parton QGP description can only be expected to apply far above that region 
as produced at colhder energies. 

At RHIC with a factor of ten increase in the initial energy density, only a continuum 
parton description of the initial conditions is tenable. The dense partonic matter so 
formed also will have much more time to evolve and produce collective signatures of its 
existence. Furthermore, a factor of ten smaller wavelength (jet) probes will finally allow 
experimentalists to resolve (i.e. see) individual quark and gluon degrees of freedom of 
that plasma. Direct observation of longitudinal work, transverse azimuthal collectivity, 
time delay, step- wise J/ip suppression in Cu + Cu, and jet quenching should, among other 
signatures, allow experimentalists to measure the properties of the partonic Yukawa phase 
at RHIC. 
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